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ABSTRACT 

A lifting surface theory has been developed 
for a helicopter rotor in forward flight for 
compressible and incompressible flow. The 
method utilizes the concept of the linearized 
acceleration potential and makes use of the 
vortex lattice procedure. Calculations demon- 
strating the application of the method are yiven 
in terms of the lift distribution on a single 
rotor, a two-bladod rotor, and a rotor with 
swept-forward and swept-back tips. In addition, 
the lift on a rotor which is vibrating in a 
pitching mode at 4/rev is given. Compressibi- 
lity effects and interference effects for a 
two-bladed rotor are discussed. 


INTRODUCTION 

Rotating lifting surfaces are an integral 
part of the propulsive unit of every aeron- 
autical and nautical vehicle, from the 
compressor and turbine blades of jet engines, 
the pumps for rocket enyines, to propeller and 
helicopter rotors. The aerodynamics of these 
rotating elements has been under extensive study 
since the advent of the airplane and with a 
combination of experimental and analytical 
approaches, succcessful designs have been 
achieved. In many cases, two-dimensional theory 
has been used, usually modified by an assumed 
spanwise distribution, and inflow velocities. 
This paper presents a compressible, lifting 
surface method for a helicopter rotor in forward 
flight within the limits of linearized theory. 

The method is based on the concept of the 
acceleration potential, originally introduced by 
Kussner (1941). The method was first applied to 
an oscillating wing in uniform translatory 
motion including effects of compressible flow by 
Runyan and Woolston (1957). The acceleration 
potential approach has now become standard for 
the determination of the unsteady aerodynamic 
forces for flutter studies of lifting surfaces 
in rectilinear motion. 

The first use of the acceleration potential 
approach for a rotating system was made in a 
paper by Hanaoka (1962) for the loading on a 
marine propeller in incompressible flow. The 
acceleration potential has been used in the past 
in studying the propeller noise problem, but in 
all of these noise propagation cases the problem 
was specialized early in the analytical develop- 
ment to the so-called far-field case usually 


with a stationary observer, whereas the lifting 
surface theory is essentially concerned with the 
details of the near-field case for a co-moviny 
observer as well as the satisfaction of certain 
of certain boundary conditions. Runyan (1973) 
utilized the acceleration potential approach to 
obtain a solution to the oscillating propeller 
in compressible flow. Oat, (1973), has derived 
a general expression for an acceleration doublet 
for any motion. Pierce and Vaidyanathan (1983) 
have treated the helicopter rotor in forward 
flight using the method cf matched asymptotic 
expansion for the incompressible case. The 
procedure developed here involves the precise 
numerical integration over the surface of the 
rotor in a time frame. The method sets forth a 
formulation of a fundamental three dimensional, 
compressible, unsteady aerodynamic theory for 
propellers and helicopter rotors. 

The next section contains a brief 
derivation of the fundamental equations, delud- 
ing a discussion of some implications of the 
equations. The third section contains a 
description of the method of solution. Finally, 
the results of some calculations for the several 
specific examples a 1 e given. 
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SYMROLS 

rotor blade area 
aerodynamic influence 
coefficients 
Fourier coefficients 
speed of sound 
chord of rotor 

thrust coefficient per blade 
(thrust/npft 2 R t 4 ) 

vector distance from doublet to 
downwash point 
absolute value of 8 

unit vector of 8 

value of sinyular integral 

kernel function 

unit vector at downwash point, 
normal to velocity vector 
unit vector at doublet point, 
normal to velocity vector 

direction cosines of n 

direction cosines of n 0 
pressure 

position vector of doublet 
from inertial frame origin 
position vector of downwash 
point from inertial frame origin 
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source or doublet strength 
rotor tip radius 
rotor root radius 
distance of downwash point 
along the span 

distance of doublet along the 
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upper limit of spanwise panel 
lower limit of spanwise panel 

distance of doublet along 
span at singular point time 
field time 

velocity of rotor system, 
parallel to x-axis, positive in 
negative x-direction 

velocity at downwash points 
♦ 

velocity component of V at 
the downwash point normal to 
the rotor leading edge 

velocity of doublet 
velocity of rotor system, 
parallel to ; axis 
downwash velocity 
distance from pitch axis to 
downwash point 
Cartesian coordinates of 
downwash point 
Cartesian coordinates of 
doublet position 
twist angle at downwash point 
twist angle at doublet 
position 

angle of axis of rotation 
relative to z-axis 
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V 0 /c 

* . i!?o 

L C 3x 

angular position of blade at 
time t 

angular position of blade at 
time t 

blade angle of attack 

blade angle relative to plane 

of rotation 

advance ratio 

air density 

time 

time at which Integrand in Eq. 

(24) becomes singular 

velocity potential 

source acceleration potential 

doublet acceleration 

potential 

azimuth angle 

rotation speed of rotor 

vibration frequency of rotor 


BASIC FORMULATION 

The formulation of the aerodynamic 
equations is based on the linearized 
acceleration potential approach. The fluid Is 
considered perfect, with no separation and the 
formulation Is based upon the assumption of 
small perturbations. The wake created by the 


lifting rotor is assumed to lie in the skewed 
helical path taken by the rotor blade. One 
reason for adopting the acceleration potential 
approach is that the pressure discontinuity 
occurs only on the surface of the blade and thus 
the boundary conditions need only be applied on 
the blade surface and not throughout the wake. 
The blade is treated as a very thin surface of 
discontinuity across which a pressure jump 
occurs. Th effect of compressibiity is taken 
into account by utilizing the complete 
linearized potential for a lifting doublet, 
along with the effects of retarded time. 

As shown in Fig. 1, an inertial coordinate 
system has been used in which the origin of 
coordinates is fixed to a point on the yround. 
The helicopter rotor is moving in the negative 
x-direction with velocity U, in the positive 
z-direction with velocity W and is rotating 
:ounter clockwise with a constant angular 
velocity fl. A point of interest on the rotor 

+ 

blade is designated by the radius vector X q (t) 

from the origin of the ground based coordinate 
system. 

Let ¥ be the acceleration potential of a 
source (or doublet), the perturbation pressure 
is then given by 


P = -P* (1) 

This expression represents the pressure p at 
+ 

point X due to a single source 

4 

(or doublet) located at X 0 . The potential ¥ 
contains a constant "q" which represents the 
strength of the source and thus the mar^itude of 
the pressure. In this form, there Is no 
boundary condition available to deten-.ina the 
value of 'V and the resulting pressure. 

Recourse can De made to the velocity potential, 
since the spatial derivative of a velocity 
potential represents a velocity. The 
relationship between the pressure and velocity 
potential for an inertial coordinate system is 

P-*$ (2) 

where ^ Is the substantial derivative. 

Dropping out the second order terms 

and integrating with respect to field time 

results in 

t 


♦(t) - J T(t') df (3) 


The acceleration potential ¥ s satisfies 
the wave equation 

^ ~r ■ - 4 * ( 4 ) 

c at 

where f(X,t) Is a source distribution. Further- 
more, if the path of an Isolated source is a 

function of time variable, X 0 (t), then 

f(X,t) ■ <s(X - X 0 ) where* 6 Is the delta function. 
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Using the Green's function formulation the 
acceleration potential expression for a moviny 
source, f s can be written as (Morse and 
Feshbach, 1978, p. 841) 
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(b) 

4- 

where X 0 (t) designates the position of the 


source at time t, X is the position of the 

4- 

field point at the time t, V 0 (r) is the 
velocity of the source point at time t, c is 
the speed of sound and q is the strength of 
the source. An auxiliary equation which relates 
the time interval (t - t) to the distance 
between the two points is 


t 


| * - *p(T) 

c 


( 6 ) 


which is usually referred to as the causality 
condition. Eq. (5) expresse: the potential as 
an explicit function of t, and only through Eq. 

(6) as an implicit function of t ans t. From 
Eq. (3), the velocity potential due to a moving 
source is 


♦ s (t) 


t t 

/ f (t 1 ) dt' = / 
-«0 — 00 


dt', 

4ir[n - 5 • 5j t , ( 1 


required. This derivative Is taken normal to 
the flight path at the location of the downwash 
point, as follows 


Aw = —— 5 - = i* • V* d>. 
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o 

Eq. (9) gives the downwash at a field point 
(x, y, z, t) due to a doublet placed at a point 
(*o*yo>Zo» T ) having a strength q. In 
order to represent a lifting surface such as a 
rotor, it is necessary to distribute the 
doublets over the lifting surface and integrate 
over the surface to obtain the downwash at a 
field point. If the downwash is known, the 
quantity "q" can be determined. Letting K be 
the expression on the RHS of Eq. (9), the final 
equation is 

w n - // K dA' (10) 



dT\ 


where ft * ? - jt Q , D = Iftl 

and dt' = [1 - M] dT > 

The quantities t' , t' and t, r satisfy Eq. 


By definition, the doublet velocity 

potential ^ of a doublet aligned along n 0 
can be written as 


♦n (t) = ar£ *s (t) = "o • 7 ? 0 *s = - *0 ■ 7 x 
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Note that for Incompressible flow, c the 
first term 0 and the Integral remains 
unchanged except for the upper limit where x ■ 


To obtain the final equation for downwash 
Aw n , a second directional derivative is 


where A' is the area of the rotor surface. 

The LHS, w n , represents the known boundary 
condition and is the velocity normal to the 
velocity vector at the downwash point. By using 
the no flow condition for the velocity 
perpendicular to the blade surface, the velocity 

component in the n direction is V n tan9 w or 

w n (r,t) ■ v n tane w = // K dA ‘ UD 

where V n is the velocity component of V at the 
downwash point and is normal to the rotor 
leading edge and % is the angle of attack. 

Thus the problem requires setting up a method of 
solution of Eq. (11) from which a value of q, 
the unknown doublet strength, can be determined 
which satisfy the known velocity boundary 
conditions w n . 

This represents a rather formidable 
computing task and the history of lifting 
surface theory even for non-rotating wings has 
centered on devising approximate methods to 
accomplish the integration in an economical 
manner. One method, termed the vortex lattice 
method, has been very successfully applied to 
aircraft wings, and is probably the more 
economical procedure of the many variants. This 
method was first demonstrated for the unsteady 
case by Runyan and Woolston (1957) and was later 
expanded by Albano and Rodden (1969). This is 
the method adopted in this paper and the 
application will be discussed later. 
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Specification of Coordinate System 

The blade has the cord C and length R^ - 
R s , Re being distance to the the root of 
the blade, Rt is the distance to the tip of 
the blade. Let the blade momentarily coincide 
with the coordinate system along the positive 
x-axis at t s 0 and execute a counterclockwise 
rotation with angular velocity n while moviny 
with velocity U along the negative x direction 
and velocity W along the positive z 

direction. Since the vortex lattice method has 
been adopted, the doublet point lies C/4 ahead 
and the downwash point lies C/4 aft of the 
section midchord. The position of the doublet 
point as well as the downwash point can be 
established as follows. The Cartesian 
components of the doublet position are 

x 0 = -Ut + r 0 cos(nx) - (C/4)sin(nx)cos oq 
y 0 = r 0 si n(flx)+ (C/4) cos(8x)cos a 0 
z 0 = Wx + (C/4) sin oq (12) 

where r 0 is the radial distance of the doublet 
along the span. With the substitution of 
C ♦ -C, r 0 + r, x ♦ t the position of the 
downwash point is given by 

x = -Ut + r cos(nt)+ (C/4) sin(nt)cos o 
y * r si n( fit) - (C/4) cos(nt) cos a (13) 

z = Wt - C/4 sin a 

In Eqs. (12) and (13), the angles a,ao are the 

twist angles of the velocity vectors ? and tf 0 , 
respectively, defined by 


W(U + r Q n $1n(flx)+ (C/4 )fl cos(nx)cos o Q ) 


V n / W fc + V ‘ 

O 0 


W(r Q n cos (fix) - (C/4)n sin(ftx)cos a Q ) 


m o ‘ " 


V„ /w 1 - + V / 

O 0 


V 

n = 

/I — 

/ tr + v ‘ 


V Q ' 2 - (U + r Q fl sin(8x)+ tC/4 )fi cos(nx)cos a o ) 2 
+ (r Q fi cos(fix) - ( C/4) n sin(fix)cos o Q ) 2 (17) 


By the same procedure, n = li + mj + nk, 
where 

i = W(U + rn sin(Qt)- (C/4 ) n cos(nt)cos a) 

vV W 2 + V 2 


-U(rn cos(nt)+ (C/4)n sin(nt)cos a) 


V / W c + V 


tan a = U sTn(nt)+ rn 


' an a o " U si n(fix)+ r Q n 




The reference piupp defined by the doublets and 
downwash points is a twisted surface. From 

Eq. (12) the doublet velocity can be computed, 
namely the time derivative of the position 
vectors. 


V + v + z o k 


The unit vector n 0 is chosen to be 
perpendicular to the twisted surface created 


V‘ 2 = (U + rfi si n( fit) - ( C/4) n cos(nt)cos a) 2 
+ (rn cos(nt)+ (c/4)n sin(nt)cos a) 2 (19) 


the vector D * X-X 0 defined in Eq. (7) can be 
expressed as 

D » j[u( t- t) + r cos(nt)- r Q cos(nx) 

+ (C/4)(sin(nt)cos a+ sin(nx)cos a Q )] 2 


by the velocity vector V 0 which is a function 
of r 0 , throuyh Eq. (14). 


Express n 0 as 


+ + ♦ 

3 JU + mj + n k 
oo o o 


where t f m , n are the directional cosines of 
o o o 


+ [r sin(nt)- r Q sin(nx) v 

P 

-(C/4)(cos(nt)cos o+ cos(nx)cos a Q )l’ 

+ [W(t-r) -(C/4) (sin a Q + sin a )] 2 } 1/2 

With the substitution of the quantities, the 
integral Eq. (11) was solved for the unknown 
q(r 0 ,x) by using a collocation process based 


the unit vector ng . It can be shown that 


!<& 


XX 



on the vortex lattice assumption. The kernel Is 
singular when D * 0 t and this was handled by use 
of the finite part technique. 


SOLUTION OF INTEGRAL EQUATION 

In following the vortex lattice technique 
the rotor Is divided Into a number of 
predetermined panels, both spanwlse and 
chordwlse, In each chordwlse panel, a line of 
doublets of unknown strength q-j Is located at 
the 25% chordwlse location of the particular 
panel, and the downwash is evaluated at the 
point located at 75% chordwlse location of the 
panel. Therefore, a collocation procedure is 
used to obtain a set of equations in terms of 
the unknown loadings qi . It is also assumed 
that the spanwlse loading q-j Is constant along 
each of the panels. A set of equations is thus 
obtained as shown below. 

w n ■ I A nm q m < 21 > 

m 

r 

where Afim * / u dr and where n refers to 
nm ' nm o 


singularity. The integration domain was divided 
Into areas as shown in Fig. 2. Areas 1-4 
(hatched) were computed numerically using a 
two-dimensional Romberg integration (Davis and 
Pabinowitz, 1967) and the contribution of the 
singular region {unhatched) was obtained in 
closed form by consideration of the finite part 
as shown in the next section. 

Treatment of Singular Term in Integral - The 
integral in the downwash equation, Eq. (11), is 
singular when D+O and produces a complication 
which must be treated properly. It should be 
remembered that the integration path along "x" 

Is the path the doublet has taken in arriving at 
the final doublet point at (c/4, r 0 ) measured 
in the local blade coordinates and can be 
considered as the wake. The Integration takes 
place along the path from -• to the final 
doublet position at t 0 . The distance D is the 
distance from the integration point at time t to 

the downwash point at £• 

There is a particular set of values of r 0 
and t for which the denominator D approaches 
zero, thus resulting in an infinite integrand. 
The singular part of the Eq. (11) is 


the downwash point and m refers to the vortex 
lattice. The kernel K is a complicated function 
which Involves an integration over t. 


j u i 2 

A T 1 


n-n„ - 3(0 *n) (D • nj 


dx dr 0 (24) 


The term q(r 0 ,x) represents the strength 
of the doublet located at r 0 and at time x, 
and Is proportional to the unknown loading. In 
order to account for unsteadiness, a solution 
was formulated to take into account the time 
variation of the strength of the wake. This was 
done by assuming a Fourier series of the form 
m 

q(r Q ,T) « A q + l (A n co$(nftx) + B n sin(nftx}) (22) 

If q(r 0 ,x) Is assumed to be a function of r 0 
alone, which means that the wake strength does 
not vary with time, the Fourier series reduces 
to q(r 0 ) * Aq. A solution obtained with 
this approximation Is termed the quasi - 
steady solution. 

This series was Inserted in the basic 
equation and Integrated with respect to t. 
However, there were more unknowns than 
simultaneous equations to solve for the 
unknowns. The additional required equations 
were obtained by evaluating Eq. (11) at a number 
of azimuth locations. For Instance if m ■ 1, 
then 

q(r 0# t) * *o + A 1 C0St,T + B i s1n0x * (23) 

The azimuth was divided Into equal segments of 
120° and the proper boundary conditions 
applied at f ■ 0°, 120°, and 240° thus 
providing the necessary additional equations. 

Numerical Integration of Kernel 

The Integration was performed hy numerical 
Integration, except for the area surrounding the 


As 5*0 at the downwash point, D becomes 

perpendicular to n, therefore, at the singular 
point, the second term is zero and will be 
neglected in the treatment of the singularity. 
However, this second term is retained in all of 
the numerical integrations involving Areas 1-4 
since it represents an important contribution 
particularly when the blade Is passing over a 
trailing wake. 

The time and distance at which the integral I 

A A 

becomes singular are designated by x and r 0 . 

The domain of the integration in Eq. (24) 
consists of a rectangle in which the duration 
x 2 -xi is kept extremely small. In other 
words, the integration Is performed along a slit 
in r 0 , over which the 2nd term in Eq. (24) is 
negligible. Therefore the Integral I can be 
approximated by 

r M x 9 n • n rt 

I ■ IS 2 — g — - di dr (25) 
r t T 1 D 

Furthermore, noticing that D 2 is quadratic in 
r 0 , if on is independent of r 0 , then the 
Integration on r 0 can be performed 
analytically. This can be achieved by 
recognizing that in the vortex lattice method, 
the rotor is divided into spanwlse panels from 
r A to r u . If these spanwlse panels are 
small then the variation In oq Is small. 

da o WO 

*0 (U Slne 0 + r Q o) Z + W 2 * 26 ^ 

If the value of oq Is approximated by Its mid 
panel value, It Is possible to Integrate Eq. 
(24), In closed form In the r 0 direction. 

This Is quite acceptable In the helicopter mode, 
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because doo/dr 0 is in the order of 

magnitude 10“3 or smaller. The value n 0 is 
also a function of r 0 , hut in the region of 
the singularity it has a very small variation 
and is evaluated at the singular position, 

(T,r 0 ). Performing the r 0 Integration 
results in the form 



where <j(x) is a function containing all the 
non-singular part after performing the r 0 
integration and t(x) * 0, at 

t=t (ti < t < 13). It: can be argued 
physically that since the quantity D(x,r 0 ; 
t,r) as well as its modified form f ( x) (after 
tegration over r 0 ) represents the distance 
between two points in space it must be positive 
and real for all its arguments, and ne/er become 
negative. Denote the value of r 0 and t at 

A A 

which D becomes zero as r 0 and t 0 . Thus, in 

the neighborhood of t the function f(*r) behaves 
like a parabolic function and has a second order 
zero. 

Expanding f( t) in a Taylor series about the 
singular point t results in 

f ( t ) = f{T) + + f'(0 (t - Wn +... 

(28) 

Since t is a second order zero 

and f(i) = f'(i) - 0 (29) 

Eq. (29) has been verified numerically. If only 
the square term is kept in Eq. (28), Eq. (27) 
can be written as 

i * l} -4- 1-44- + aL4L ♦ dT < 30 > 

T 1 f“(t) (t-t ) i 2 (t-t) 2 


as At Is kept large because the very large 
values of the integrand near the singularity are 
avoided. On the other hand, regarding the 
finite part integration, the denominator was 
expanded in a Taylor series about the 

singular point, x. Therefore, it is desirable 
to maintain Ax as small as*possible to keep 
within the limits of the applicability of the 
series expansion. Numerous calculations were 
made, varying Ax until a reasonable 
convergence was found. This value was found to 

be .01 (t- t) , i.e. 1% of, the time difference. 
Actually, there is very little difference 
between 1% or 10% of the time difference and the 
computing time and cost is considerably reduced 
by using 10%. For trend studies 10% is 
recommended principally to reduce computer 
costs. However, for final design type analysis, 
a smaller value of time difference Ax is more 
appropriate. 

For the spanwise direction, Ar 0 is also 
an integration limit variable. The finite part 
integral was obtained by approximating the angle 
of twist of the velocity vector across a segment 
by assuming it constant across the segment, 
having a value as determined at the center of 
segment. Numerical experimentation indicates 
that for a helicopter, Ar 0 = 0 is satis- 
factory. 


APPLICATION TO SPECIFIC EXAMPLES 

The foregoing analysis has been applied to 
several specific examples which are given in 
Figs. (3) and (4). The following section 
presents results for several paneling 
configurations; e.g. 5 spanwise and 1 chordwise 
panels (designated (5-1)) and 7 spanwise and 3 
chordwise (designated (7-3)). The rotor blade 
was maintained at a constant pitch setting of 
% * .1 radians for all the calculations. 


In Eq. (30), if x2 and are chosen 


Single Blade 


symmetrically about x, then the odd derivative 
terms Integrate to zero. Futhermore, the third 

A 

te^m can be neglected since g M (t) is small. The 
major contribution comes from the first term. 
Then using the standard integration technique 
(Mangier, ' 1952) the final result for the 
integral is 


i . . 

f"(r) 


(31) 


A 

where 2Ax ■ X2 - and ti < t < , 

A numerical problem arises because the 
finite part integration results in a negative 
number which is close to the total of the 
surrounding numerical integration areas which 
are positive. Thus, It Is necessary to take the 
difference between large numbers, and the final 
Integration accuracy is dependent on the 
accuracy of the two Integrations. On the one 
hand, the numerical integration is more accurate 


In order to investigate the convergence of 
the method when using the vortex lattice 
procedure, the program was run for several 
chordwise and spanwise elements for the incom- 
pressible case. The thrust coefficient Ct 
vs. the azimuth angle is shown in fly, (5). (In 
all of the following plots for thrust coeffi- 
cient vs. azimuth angle, the thrust was 
calculated for 16 uniformly spaced azimuth 
angles and each curve was faired using a 
cubic spline). The rotor was first divided into 
5 spanwise and one chordwise (L*l) panel and the 
results are shown by the solid line. The 
chordwise division was increased to (5-2) and 
the results are shown by the long dashed line. 

It can be seen that very little change has taken 
place* The spanwise divisions were increased to 
(7-1) and the largest change occurred at 
* ■ 0° where the difference in Ct is about 
11%. Increasing the chordwise divisions to 3 
(7-3) shows convergence of the (7-1) case to be 
very good. 

An 1nt?rcct1ny phenomena occurs In the 
reylon of small azimuth angles. For \jmO to 
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37°, the lift Increases to a local maximum at 
♦*370 then the lift abruptly falls to a local 
minimum for *»60° and then rapidly Increased 
to a maximum at A similar phenomenon 

Is shown analytically by Eglof and Landgrebe 
(1983) In Fig. 60 of that report where a 
local minimum and a local maximum occur In the 
same range of azimuth angles, even though the 
geometry of the two blades and the flight 
conditions are different. Also, In Fig. 93 of 
the same report some test data shows a similar 
variation of loading In the same azimuth range. 

The chordwlse pressure distributions for 
the (7-3) case are presented in figure 6. It 
should be remembered that in using the vortex 
lattice method, the loading is concentrated at 
the location of the vortex which for the (7-3) 
case Is located at .0833C, .4160, and .750. The 
pressure was faired using a cubic spline through 
the three vortex locations and the known value 
of zero at the trailing edge. The distributions 
are given for 7 spanwise positions. In general, 
the curves exhibit the expected shape, having 
the largest values as the leading edge is 
approached. For the span distribution the 
values at r/Rj - .85 are slightly larger than 
the values at r/RT * .95,, Indicating a falling 
off In the tip region. 

From these concentrated forces, the section 
pitching moment can be calculated. Figure 7 
presents these results for ♦ * 90 degrees. The 
section moment was taken about the 1/4 C and a 
nose down moment is taken as positive. The 
pitching moment shows some rather dramatic 
changes along the span. The moment is nose up 
near the tip ( r/R j * .95), changes to a small 
nose down value, then becomes nose up for most 
of the Inboard region. Integration of the 
moment would result in a total pitch moment up 
at * * 90°. 


Swept Tip 

The segments used for the vortex lattice 
for the swept tip studies were(5-l), where two 
equal segments were used In the tip region and 
three equal segments were used In the unswept 
Inboard section. In Fig. 8 the lift Is shown 
plotted against azimuth for the two sweep 
conditions and for zero sweep. In general, the 
three results show little difference. The 
sweptback configuration has a larger lift from 
♦ • 300° to 400. For ♦ * IQOO to 2400, 
the swept forward configuration has a slfyntly 
larger lift. It appears that the total lift for 
one rotation for the swept-back case and the 
swept forward case would give about the same lift 
as produced by the unswept rotor. In Fig. 9 
the lift distribution along the totor span Is 
given for 0°. The major effect of sweep 
Is concentrated at the tip, where the swept-back 
tip load 1$ greater than both the unswept and 
sweptback cases . In fig, lo, 180°. 
Comparing to fly. 9, the swept-back tip load is 
larger than both the unswept and the 
swept-forward tips. 


Blade Oscillating in Pitch 

An example of unsteady loads on a rotor 
blade with (5-1) paneling which Is oscillating 
In a pitching mode about the mid-chord at a 
frequency of 4 per revolution (120 cycles/sec) 
is yiven on fig. 11. For this case a 17 term 
Fourier series (m=8) was used to simulate the 
oscillating load, which was comprised of one 
constant term, 8 cosine terms, and 8 sine 
terms. The steady and unsteady rotor blade 
loading is yiven for one revolution. The blade 
was oscillated through an angle of .1 rad. 
about a mean anyle of .1 rad. The effect of the 
oscillation Is readily apparent as compared to 
;he steady case. With the harmonic 
representation of the loading, the magnitude and 
phase of the several harmonic loads are easily 
determined. The magnitudes are plotted in 
Fig. 12. The only harmonic loads that were 
significantly changed from the steady case were 
the 3rd, 4th and 5th. Both the 3rd and 5th 
harmonics were Inc- **d and the 4th harmonic 
was dramatically incr.'sed. Another calculation 
was made for the non-osci 1 latory unsteady case 
and compared to the quasi -steady case. 

Virtually no difference was observed, indicating 
that, at least for this case, the rate of change 
of loading in a revolution of the blade is small 
enough so that the effect of a variable wake is 
negligible. 


Compressible Effects (5-1) 

For a one-bladed rotor, the effect of 
compressibility is illustrated in Fig. 13, in 
which the Cj is plotted against azimuth 
angle. The Incompressible result is Included 
for comparison. As expected, the compressible 
load is larger than the incompressible 
throughout one revolution. The effect is 
greatest In the reyion of the advancing blade 
and smallest in the retreating region as would 
be expected. 


Two-Bladed Rotor In Compressible Flow (5-1 per 
blade) 

The method has been extended to the 
two-bladed rotor for the compressible case and 
the results are shown in Fig. 14. The thrust 
coefficient Cj per olada is given vs. azimuth 
angle for a single bladed rotor and for a 
two-bladed rotor. For azimuth angles from 
t ■ 20° to 120° the single blade rotor has a 
larger Ct* For ♦ • 120° to 260°, the Cj 
on the one and two-bladed rotors are 
approximately the same. However, for 260° 
to 340° a dramatic reduction in lift occurs 
for the two-bladed rotor as compared to the one 
bladed results. The lowest lift occurs at 
ik* 292° which places the other blade of the 
two-bladed rotor at ♦ * 112°, the point of 
maximum lift on the other blade. Apparently the 
high lift on the blade at 112° creates a 
very unfavorable Induced velocity on the second 
blade at i(j* 292° which requires the loading 
to go to zero In order to satisfy the boundary 
conditions at 292°. 
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CONCLUDING REMARKS 

A linearized lifting surface theory 
including the effects of compressibility has 
been developed for a helicopter rotor in forward 
flight. The method utilizes the concept of the 
acceleration potential, and makes use of the 
vortex-lattice procedure for performing the 
required Integrations. In addition, the method 
has been extended to include the effects of 
unsteady flow. 

Sample calculations have been done for 
several cases. These include the effect of 
swept-back and swept- forward tip. The effect of 
these two tip configurations was minimal on the 
total loading for one revolution. However, the 
loading distribution changed considerably for 
several azimuth positions. A comparison of the 
thrust coefficient, Cf, of a one bladed rotor 
and a two bladed rotor was made. In the 
azimuthal range between 20° and 120°, the 
one bladed rotor showed higher lift. However 
between ^ s 260° to 340° the two bladed 
rotor indicated a lower Cf. Compressibility 
was Investigated for one configuration. As 
expected, the effect was greatest in the 
advancing blade region ( 900) and was 

minimal in the retreating blade region. The 
effect on Cj of a bladu oscillating In pitch 
at 4/ rev Is given. The effect on the total 
blade lift Is shown and the effect of the 
oscillation is readily apparent. The harmonic 
content was calculated and the greatest 
difference between the oscillatory and 
non-oscillatory cases was found in the 4th 
harmonic. 
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Fig. 2 Integration Areas 
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Fig. 3 Unswept Configuration and nput 
Parameters 


Fig. 4 Swept Tip Configurations 




Fig. 5 Thrust Coefficient vs. Azimuth Angle for 
Four Panel Configurations for a Single 
Rotor Blade, Incompressible, (u ■ 0.17, 
Ojj * 0.1 rad, eg. » .OF rad, 

O ■ 30 rad/sec) 


Fig. 6 Chordwlse Pressure Distribution for 
Several Spanwlse Locations, p » 90°, 
incompressible, (p • .17, e 9 « 0.1 
rad, og. » .05 rad, a ■ 30 rad/sec) 



Fig. 7 Spanwlse Section Moment Distribution 
about 1/4C - Positive Nnse Down, 

■ 90°, incompressible, 

(u ■ 0.17, Og ■ 0.1 rad, eg. ■ .05 
rad, Q ■ 30 rad/ sec) 


Fig. 8 Comparison of Lift on a Swept-Back Zero 
Sweep and Swept-Forward Blade, 
Incompressible, (u a o.l7, e. - o.l 
rad, op a ,05 rad, o ■ 30 rad/ sec) 
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Fig. 9 Spanwlse Section Lift Distribution for Fig. 10 Spanwlse Section Lift Distribution for 

Swept-TIp Configurations, * » O'*, Swept-TIp Configurations, * ■ 1800, 

incompressible, (u ■ 0.17, 0« « 0.1 Incompressible, (u • 0.17, e« ■ 0.1 

rad, a,- » .05 rad, n » 30 rad/sec) rad, op ■ .05 rad, Q ■ 30 rad/sec) 
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Fig. 11 Comparison of Lift on a Rotor Blado 
Oscillating In Pitch at 4/Rev. to tha 
Lift on a Kon-Oselllatory blado, 
Incomprosslblo (it • 0.17, • 0.1 

rad, op * .05 rad, Q ■ 3U rad/ sec) 


Fig. 12 Harmonic Content for Non-Osclllatory and 
Oscillator/ Cases - r/;»r ■ .95, 
incompressible, (it ■ 0.17, e* ■ 0.1 
rad, op * <05 rad, 0 • 30 rad, sec) 
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Fig. 13 Incompressible and Compressible Lift for 
One-Bladed Rotor, (y * 0.17, 8« = 0.1 
rad, op * .05 rad, Mjip 3 0.54} 



Angle, deg 


Fig. 14 Lift on Two-Bladed and One-Bladed Rotor 
vs. Azimuth Ang’e, compressible, 

(y * 0.17, 0g 3 0.1 rad, ap = .05 
rad, Kf jp * 0.54) 





DISCUSSION 
Paper No. 7 

LIFTING SURFACE THEORY FOR A HELICOPTER ROTOR IN FORWARD FLIGHT 

H. Tai 
and 

Harry L. Runyan 


Bob Sopher, Sikorsky Aircraft : Just a couple of questions. One, you had a rigid rotor in this 

analysis? 

Tai : Yes. 

Sopher : There is no aeroelasticity and you went up to a » of .17? 

Tai: That's correct. 

Sopher : First of all, I think that as soon as you put aeroelasticity in you are going to see 

different trends when you put sweep on in comparison without sweep. 

Tai: Oh, sure. 

Sopher : The second thing is at u equal to .17 I think that the skewed helical wake will not 
be a very good approximation to the actual wake, which will be substantially distorted— more 
like a wake that you would get under hover conditions. So I question the utility of assuming a 
skewed hi ileal wake. 

Tai : As T pointed out this is a first cut. Of course later you can modify this by iteration 

processes. You can go baok and calculate the wake and put it back again. Hopefully, that 
procedure would ^ve you a better result. 

Sopher : The third question is what advance ratios do you expect to apply the analysis up to? 

Tai: The answer is I don't know. However, I would think that the higher forward speed would 

probably have a better answer. Because you don't depend on the wake that much. 

Sopher : I question that because as you go up to higher speeds you are going to find that you 
run into situations where you get transonic flow on the advancing blade and I do not believe 
that the linear analysis will apply accurately under those conditions. 

Tai : Perhaps you are right. 

Sopher : As a matter of fact this research center has developed transonic flow analyses which 

apply to three dimensional lifting blades so I would say that the primary utility that I would 
see in this analysis is for hover applications where the linear analysis is valid, but you would 
have to use a distorted wake. 

Tai: For hover cases you would really expect the wakes to stack up and then you end up with a 
very difficult mathematical problem. However, I guess from my past experience you probably can 
get the loading by some numerical procedure. For example, you can do extrapolation. Assuming a 
certain W and then you extrapolate for W s 0. I don't know. We don't have a dear under- 
standing. I admit that. 

Jim MoCroakey. U.S. Army Aeromechanics Laboratory : You have made some nice progress on this 

approach since you talked with us a year or so ago. It's interesting; it's nloe to see some 

results being generated for some revUstic oases. I wanted to ask a couple of minor questions. 

I presume this blade is untwisted, is that right? 

Tai: Yes. 

MoCroskey : How did you treat the reverse flow region? 

Tai : The twist can be added on very easily beoause [we only have to] add on the boundary condi- 
tion. To avoid a reverse phenomenon we deliberately use very large cutoff. You can see [that 
it is]_six feet. We try to avoid that region. 

MoCroskev : In faot the sketoh in the book is a little misleading beoause the cutout is more 

like 30$ instead of the 10 or 15J [that appears in the figure]. So you Just avoided it by 

having a root out out. 

Tai: Yea. 





MoCroakev : The final thing la on the influence of compressibility. Do you have some Idea of 
what you predict in Figure 13; is that a larger effect than you would predict if you Just used 
some kind of Prandtl-Glauert scaling on the incompressible solution? 

Tai : Well, I guess we should be blamed for not making the abstract very clear. We used a very, 

very honest way of doing that. We didn't use rny approximation at all. In other words, as I 
pointed out, you find the t as function of r Q . In other words, you give it the radius of the 
doublet [and] you go there and find the t which serves as your upper limit— which is not a 
trivial matter. To answer your question we say that we use true, honest compressibility effect. 

McCroskey : But the question is how good would the Prandtl-Glauert type approximation be to what 
you actually calculated? 

Tai : Well, to honestly answer that question — we don't know. We didn't check [it], but I think 

it is not very easy to check it out. 

Bob Ormiston. U.S. Army Aeromechanics Laboratory : I want to commend your results. It looks 
like you have made some pretty good progress in the last year or so. The methods we are using 
now for routine rotor loads analysis are usually based on some fairly primitive assumptions like 
strip theory and 2-D airfoil coefficients and so forth and what we ultimately have to get to is 
very, very sophisticated, maybe, 3-D CFD kinds of analyses. It looks like what you've got is an 
intermediate type of analysis which may be very practical. Hy question is do you think there is 
a practical way to generalize the results you have gotten, say, to come up with generalized 
forcing functions for specific loading distributions, a family of loading distributions that you 
might be able to calculate and then not have to repeat the integration problem for each particu- 
lar configuration that you are analyzing? Is there a practical way to do that? 

Tai : Well to answer your question, the answer is yes. I did not mention that when we break the 
blade into different segments. Apparently the matrix is highly diagonal. In other words off 
diagonal matrix you can use less accurate methods to generate. Beyond that, to answer your 
question, I think in a practical sense we can generate those matrix elements and store them and 
only change the boundary conditions to do all the types of calculations. In other words, the 
answer is indeed it can be very practical. 
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